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Gravitational-wave data analysis requires a detailed understanding of the highly relativistic, late 
stages of inspiral of neutron-star and black-hole binaries. A promising method to compute the late 
inspiral and its emitted waves is numerical relativity in co-rotating coordinates. The coordinates 
must be kept co-rotating via an appropriate choice of numerical relativity’s lapse and shift functions. 
This article proposes a model problem for testing the ability of various lapse and shift prescriptions 
to keep the coordinates co-rotating. 
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I. INTRODUCTION 

Patrick Brady, Jolien Creighton and I Q have dis¬ 
cussed the importance, for gravitational-wave data anal¬ 
ysis, of understanding the evolution of inspiraling black- 
hole binaries. We showed that there is a crucial gap be¬ 
tween the point, during inspiral, when a post-Newtonian 
analysis breaks down, and the point, just before merger, 
when conventional numerical relativity techniques can 
take over the analysis. We argued that this IBBH (In¬ 
termediate Binary Black Hole) gap might best be filled 
by numerical relativity computations carried out using 
co-rotating coordinates. In such coordinates, the quanti¬ 
ties being computed (components of the 3-metric ^ij and 
extrinsic curvature Kij ) change on the timescale of inspi¬ 
ral T*, which is much longer than the orbital timescale 
1/Cl — 1/(orbital angular velocity). Therefore, in such 
coordinates one may be able to take long time steps, con¬ 
strained only by r*, thereby speeding up the computer 
computations by a large factor relative to those carried 
out in the usual asymptotically inertial coordinates. 

Numerical relativity in co-rotating coordinates may 
also be effective, and perhaps necessary, for analyses of 
the late stages of inspiral of neutron-star binaries and 
neutron-star / black-hole binaries. 

A key foundation for such computations is a prescrip¬ 
tion for choosing the lapse and shift functions a and /?*, 
which generate the coordinate system as the computa¬ 
tion proceeds. The shift function P’’ must keep the coor¬ 
dinates co-rotating with the binary’s orbital motion, and 
the lapse function a must appropriately retard the time 
slicing in regions of strengthening gravitational redshift; 
and this must be done quantitatively in such a way that 
the combined effects of a and /3® keep djij/dt ~ 
and dKij/dt ~ For simplicity, I refer to such co¬ 

ordinates as co-rotating, even though they require more 
than just co-rotation. 

Brady, Creighton and I have proposed a class of 
equations for the lapse and shift, any one of which (we ar¬ 
gued) may keep the coordinates co-rotating, if they have 


been chosen co-rotating to begin with. The simplest of 
these are the “minimal-strain” lapse and shift equations: 


DJ(D(,/3,)-aK,,) = 0 , (la) 

in which a is to be replaced by the following expression 
in terms of Pj 




(lb) 


Here Di is the covariant derivative in the slice of constant 
time t, which has 3-metric 7 ^ and extrinsic curvature 

Numerical relativists, in conversations with me, have 
suggested other choices of lapse and shift that might keep 
the coordinates co-rotating. A detailed tests of their 
proposals and of Brady, Creighton’s and mine are now 
needed. The purpose of this Brief Report is to suggest 
a model problem that can serve as a test bed for the 
various proposed lapse and shift equations. 

This test bed is a specific 4-dimensional spacetime that 
resembles that of a fully relativistic, inspiraling compact 
binary. More specifically, it is a spacetime with a 4- 
metric that possesses, semi-quantitatively correctly, all 
those features that are likely to significantly influence the 
behavior of the proposed lapse and shift equations—with 
one exception: the test-bed 4-metric does not possess 
black-hole horizons. Instead, it models a binary neutron- 
star system. This is because horizons may seriously com¬ 
plicate the task of implementing the proposed lapse and 
shift equations; so in initial tests it is best to omit them. 
It should not be difficult to generalize my test-bed 4- 
metric to include horizons, both nonrotating and rotat¬ 
ing. 

Any lapse and shift equations that stably generate 
and maintain a co-rotating coordinate system in my 
test-bed spacetime will very probably do so also in a 
real numerical-relativity calculation—i.e., in a calcula¬ 
tion that is simultaneously evolving the binary’s true 
spacetime geometry along with the coordinate system. 
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II. THE TEST-BED SPACETIME 

I shall describe the test-bed spacetime in a co-rotating 
coordinate system {T, X,Y,Z} that must not be confused 
with the coordinates {t, x, y, z} generated by some cho¬ 
sen lapse and shift equations. Thinking of the spatial 
coordinates {X,Y,Z} as Cartesian in some flat back¬ 
ground metric, we can introduce the flat-space distance 
zu from the binary’s rotation axis, spherical polar coordi¬ 
nates {R, 0, (j)), and unit vectors 9 and 0 along the 9 and 
(j) directions: 


here 2ao is the separation at time T = 0, when the re¬ 
maining time to merger is AT = t*o/ 4, and the inspiral 
timescale is r*o = ao{dao/dT)~^. 

The test-bed 4-metric is not actually a solution of Ein¬ 
stein’s field equations, but rather is cooked up to resem¬ 
ble the solutions that we expect numerical calculations 
to reveal. 

We can write the 4-metric in 3-1-1 form in the co¬ 
rotating coordinates: 

ds^ = -A^dT^ + T,j{dX^ + B^dT){dX^ + B^dT) . (7) 


= \/a2-p y2 ^ R = , 

9 = aicsm{w/R) , (j) = arctan(y/A) ; 

~ dg -wdz + Zd^ Xdx+Ydy 

fc* = — =- - - where =- 

R R ZD 

^ _ c>4, _ Xdy —Ydx 

'UJ 


( 2 ) 


Here {dx^dy^ dz} are the usual Cartesian coordinate ba¬ 
sis vectors, {dg,dfj)} are the usual spherical coordinate 
basis vectors, and is the unit vector pointing away 
from the rotation axis in the flat background metric. 

The test-bed spacetime describes two identical stars 
in a circular orbit around each other. The space¬ 
time’s 4-metric depends on four parameters: the time- 
independent mass m of each star as measured gravita¬ 
tionally far from the binary, the time-independent stel¬ 
lar radius b and slowly evolving stellar separation 2a as 
measured in the flat background metric, and the slowly 
evolving orbital angular velocity as measured far from 
the binary. 

We require that 


2m %b a , Ha ^ 1 , (3) 


so the stars are larger than their Schwarzschild radii, 
their separation is larger than the sum of their radii, and 
their orbital speeds are less than the speed of light. The 
“approximately less than” rather than strictly less than 
is because a and b are parameters measured in the flat 
background metric, not in the physical metric. We also 
require that 


da a drt 

so the inspiral is slow, 
might want to set 


— , where Hr, 1 , (4) 

r* 

To be somewhat realistic, one 


H = 



( 5 ) 


(Kepler’s law) and let the separation 2a evolve in the 
manner dictated by the quadrupolar description of radi¬ 
ation reaction Q 

a(T) = Oo (1-) , where r*o =-r—T ; (6) 

\ Tt,o J 8 


Here the co-rotating lapse A, shift and 3-metric T^ 
are not to be confused with the lapse a and shift j3 that 
are generated by some proposed prescription for keeping 
the coordinates (nearly) co-rotating, and the 3-metric 7 ij 
of that prescription. 

The co-rotating A, B® and Ty will depend explicitly 
on the spatial coordinates {A, K, Z} and will depend on 
time T solely through the slowly changing separation 2a 
and angular velocity H. 

To take account of the fact that information about 
the inspiral cannot propagate faster than light, in the 
4-metric we shall regard a and H as functions of (approx¬ 
imately) retarded time T — R rather than time T : 

a = a{T-R) , H = H(r - R) . (8) 


In numerical tests this might not be very important, since 
the the full coordinate grid is not likely to be larger than 
a few gravity-wave wavelengths, ~ 10/H, and thus will 
likely be much smaller than t* (the timescale for changes 
of a and H), except possibly near the end of a test. 

We choose the (truly) co-rotating lapse A to be 1 minus 
the (approximate) Newtonian potential of the binary, in 
accord with the Newtonian limit of general relativity: 


A = 1 — mFy , 


( 9 ) 


where the function B+ (“Newtonian potential per unit 
mass”) and an associated function B_ are defined by 


F± = 


1 


[(A-a)2-hy2 + ^2 + 52]i/2 
1 


± 


[(A + a) 2 +y 2 +^ 2 + 52 ]l /2 ’ 


( 10 ) 


and a = a(T — R). The b^ rounds off the growth of the 
potential F± as one moves inside each star. 

If our coordinates were inertial rather than co-rotating, 
the shift B would consist solely of a piece Bdrag produced 
by the dragging of inertial frames due to the binary’s or¬ 
bital motion. The transformation (j) = (f> + J fldT, from 
the inertial angular coordinate (j) to the co-rotating angu¬ 
lar coordinate (f>, augments this frame-dragging shift by 
^d(j, — flzucj), thereby giving rise to our following choice 
for the shift in co-rotating coordinates: 


2 











B = 



Ama^ 


(i?2 + a2)3/2 


Vlwcj) 


16 ma^r 2 ^ ^ 

(i?2 + a2)2-^-‘^^ ■ 

( 11 ) 


Our choice here of B^rag consists of two pieces. The 
first piece, — 4 TOa 2 (i ?2 _|_ dominates at large 

radii R ^ a, where it has the standard form and value 
— (2J/R^)d^ for the dragging of inertial frames by the 
binary’s total angular momentum J = 2ma^n. The sec¬ 
ond piece, —16ma^fl{R^ + a?)~^F-dy dominates in the 
vicinity of each star, where it has the standard form and 
value —4p/(distance to star) for the frame-dragging field 
produced by the star’s linear momentum p = Azma^dy 
(-1- for the star at x = -|-a; — for the star at —a). Because 
the vector fields dy and = Xdy — Ydx are both 

regular on the rotation axis and at the origin, our chosen 
shift B is regular. 

We split our co-rotating 3-metric into two parts: 

r., =rg + rTT. (12) 


Here 


= + (13) 

is a conformally flat part with the form (flat metric) x (1 
- Newtonian potential)2 expected from the Newtonian 
limit of general relativity; 


pTT 


-Ama^n^inRY 

(1 -hH2i?2)5/2 


(1 -I- cos^ 0) cos{2(j) -\- 2VtR){69 — ^(j)) 


-\-2 cos 9 sin(2(/) -|- 2Sli?) {9^ + ^9) 


(14) 


it represents outgoing waves rather than ingoing shows 
up in the direction of the spiral and equivalently in the 
relation {da — fld^)hy ~ hyjr^ <C flh+ and similarly for 
hx ■ An ingoing-wave pattern would exhibit the opposite 
direction of spiral and would have (dR+fldd,)hy ^ hy/r^,. 

The specific functional form of r"'’"’’, Eq. (|l^, is de¬ 
signed to have several properties: (i) It produces the de¬ 
sired radiation field ( [l^ at i? ^ 1/H. (ii) It becomes 
very small in the near zone, R ^ 1/H—negligibly small 
compared to the conformally flat part of the 3-metric 
r*^. (iii) It is regular on the rotation axis, because as one 
approaches the axis cos(2(^ — 2(j)o){99 — (j)(j>) ± sin(2(() — 
2(j)o){9<j) + (j)9) becomes a simple tidal stretch along the 
4>o direction, (iv) It is regular at the origin, because 
its angular dependence is that of the electric-type ten¬ 
sor spherical harmonic of order I = 2, which 

is constructed via Clebsch-Gordon coupling from scalar 
harmonics of order V = 0,2 and 4 Q, and the R'^ depen¬ 
dence at the origin, multiplied by these scalar harmonics, 
produces regularity. 

There is no attempt, in the chosen expression for 
to represent, even approximately correctly, the 
radiation-reaction field inside the binary. This is because 
the radiation-reaction field is so small inside the binary, 
^ m/R ||], that it is unlikely to have any 

significant influence on the behaviors of proposed lapses 
and shifts. Similarly, there is no attempt to approximate 
the details of the transition from near-zone field to ra¬ 
diation field. Except near the endpoint of inspiral, that 
transition field is small compared to m/R, the important 
part of r®. 

In summary, our test-bed 4-metric (^ has the co¬ 
rotating lapse (^), shift ([Tl|), and 3-metric (g)-(0. 


is a transverse-traceless (TT) part that represents out¬ 
going gravitational waves in the radiation zone; and two 
vectors written side by side, e.g. 9(j), represent a tensor 
product. 

The TT part of the 3-metric, requires discussion: 
In the radiation zone, R ^ A/VL, after transformation 
from the co-rotating angular coordinate (j) to the inertial- 
frame angular coordinate (j) = (/ + J ^dT ~ (f) + flT, 
takes the standard quadrupole-moment form for outgo¬ 
ing gravitational waves produced by a nearly Newtonian 
binary: F"'""'" = h+{99 — (j)^ + hx{9(j) + (j)9). Here the 
waveforms and hx are [0 

_4?77/7^0^ ^ — 

/i+ =---(1 + cos^ 9) cos[20 — 2il(T — i?)] , 

R 

—— 

hx = --2cos6»sin[2<()-2H(T-i?)] . (15) 

R 

In co-rotating coordinates (p = (f) — flT, the time oscilla¬ 
tion of this radiation field is gone. Instead of oscillating 
in time, it consists of a simple spiral pattern in space that 
changes on the slow inspiral timescale t*. The fact that 


III. TESTING PROPOSED LAPSES AND SHIFTS 


To test a proposed prescription for choosing the lapse a 
and shift —for example, the Brady-Creighton-Thorne 
minimal-strain prescription (Eqs. (|^) and associated 
boundary conditions discussed in Ref. [|I|)— one can pro¬ 
ceed as follows: 

Choose the 3-surface T = 0 as an initial 3-slice in 
spacetime, and label it with the time coordinate t = 0. 
On this initial slice, set the spatial coordinates {x,y,z} 
equal to the co-rotating coordinates {X,Y,Z}. Corre¬ 
spondingly, the initial 3-metric and extrinsic curvature 
in the {x, y, z} coordinates will be 

(i = 0) = (T = 0) , (16a) 


K,,{t = 0) 


1 



2D(^iBj^ 


T=0 


(16b) 


where the T-derivative acts solely through the depen¬ 
dence of Tij on a{T — R) and H(r — R). 
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Choose the initial lapse a{t = 0) and shift = 0) to 
be as close as the proposed lapse and shift prescriptions 
allow to 


a{t = 0) = A{T = 0), j3^ = B\T = 0). (17) 

If the prescription does not allow these choices up to a 
fractional difference of order l/(nT*) <C 1, then the pre¬ 
scription is suspect. The minimal-strain prescription (|^) 
does have this property, as do all the other prescriptions 
discussed by Brady, Creighton and Thorne | ^. 

Beginning with these initial data, use the proposed 
prescription to evolve forward in time the lapse a and 
shift /?*, and the associated time slices t{T, X,Y, Z) = 
constant and the spatial coordinates {x,y,z} on those 
time slices. At each step of the evolution, compute 
the 3-metric jij and extrinsic curvature Kij of the slice 
t = constant from the test-bed 4-metric {A, S®, T^}. An 
acceptable lapse-shift prescription will keep a, (3, jij, and 
Kij all evolving on the slow timescale r*; and most likely 
will achieve this by maintaining t ~ T, a ~ A, /3® ~ B®, 
and 7 y ~ aside from fractional differences of order 
l/(nr*). 
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